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| Compact and well defined formulae for the shadow of the virtual 

photon interacting with a large nucleus at small XBj are given in the 
QCD dipole picture. Two classes of contributions are considered: (a) 
^\ • quasi-elastic interaction of the qq dipole and (b) multi-pomeron cou- 

*7h ! pling- 

a: 

Qh! 1 Introduction 

<D 
43 

We have now at our disposal the QCD dipole picture of virtual photon inter- 
actions |, H with hadr ons ||, f| which is able to account for existing data 
^ ' on proton structure function F2 (i.e. on total virtual photon-proton cross- 

section) at small XBj- It is therefore interesting to investigate whether the 
model can be also applied to nuclear targets, in particular if it can describe 
the shadowing of the virtual photons. This becomes especially interesting in 
view of a possible extention of the HERA experimental program to nuclear 
targets ||. 

In our recent paper Q we have calculated in this model the double scat- 
tering contribution to the virtual photon shadowing^. The multiple scatter- 
ing corrections, which are more difficult to handle, were merely estimated, 
however. As the resulting shadowing at small XBj turned out to be rather 

*To be published in Acta Physica Polonica B, no. 3 (1998). 

1 The subject of photon shadowing in nuclei has a vast literature. For recent reviews, 
see e.g. @ and §]. 
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substantial and not too far from experimental data, it seems worthwhile 
to pursue the problem further and to evaluate more precisely the multiple 
scattering terms. This is the purpose of the present paper. As in m (c.f. 
also H |l0|) we are dealing with two types of contributions: 

(a) The "quasi-elastic" interaction of the qq dipole (first introduced in 
the dipole picture in |jTo|| ) contributing mostly at small mass excitations of 
the virtual photon. 

(b) The "direct" interaction of the dipole cascade (corresponding to the 
so-called multiple pomeron coupling in the Regge language) which dominates 
the high-mass excitations. 

The amplitudes describing these processes were derived and discussed in 
jD, ^, ||, ||, [[(], 11]. When applied to the interaction with nuclear targets, 



however, they must be modified in two important points: 

(i) a phase difference related to the change in the longitudinal momen- 
tum of system travelling through the nucleus must be taken into account. 
Consequently, at some point it is nessesary to transform the amplitudes into 
the momentum representation. This makes the calculation substantially 
more complicated. 

(ii) Since the transverse size of the gluon cascade is much smaller than 
that of a large nucleus, we can use the standard approximation and integrate 
over the impact-parameter dependence of the elementary amplitudes, i.e. we 
need to deal with forward scattering amplitudes only. Indeed, the convolu- 
tion of the onium and the nucleus transverse profiles is approximately 

J T A (s) T(b - s) d 2 s « T A (b) J T (b') d 2 b' = f(k = 0) T A (b) , (1) 

where T{k = 0) is the forward scattering amplitude and T A (b) is the nucleus 
profile. 

This simplifies the calculation in an essential way. Our ultimate goal is 
to give an explicit expression for the shadow, defined as 

C(A,, Bi ,Q 2 ) = ^-l, (2) 
Aa tot (n) 

where at ot (A) is the total photon-nucleus cross-section and atot(n) is the 
total photon- nucleon cross-section. atot(n) has been measured Jl^] and suc- 
cessfully computed in the dipole picture Q. This will be our input. 

To compute £j t ot(^4) we have to know the forward amplitudes for (a) 
quasi-elastic scattering of the dipole and (b) direct scattering of the dipole 
cascade from a given number, N, of nucleons. We then dress them up with 



2 



the light-cone photon wave functions |l3| and with single nucleon densities 
characterizing the nuclear shape. Finally we add them up for 2 < N < A. 
This calculation is described in the next two sections. In Sec. 4 we discuss 
the obtained results and formulate our conclusions. 

2 Forward 7*-nucleus amplitude from the direct 
interaction 
of the dipole cascade 

We start by giving the final result for the forward amplitude for scattering 
on N nucleons of the dipole cascade emerging from an onium. It reads 

F^(r,r ;Y,y) = (2N - 3)!!™ V 2 VV (zi ~ ZN)mxp 

*n / ^w^w ^(;) sffy. (3) 

\Cfc-iOO / 

where 

A(X) = ^X(A), x(A) = 2^(1) - V (l - |) - V-(|) , (4) 

a is the strong coupling constant, N c is the number of colours, tq is a 
parameter of the order of the effective size of the nucleon, r is the transverse 
size of the incident onium. The remaining symbols denote: 

Y = -log'-f, » = -log^, x F ='-f, /S=ipgj„ (5) 

and M is the mass of the dipole cascade; c is a constant giving the energy 
scale, n e g is an effective number of onia in one nucleon (estimated in 
c.f. also §). 

- / \ ■/ _ \ Q 2 +M 2 

The origin of the phase factor e 1 ^ 1 z NJ mx P = e H z i z nj 2v where z\ 
and zn are the longitunal positions of the first and the last scattering 
(m is nucleon mass and v is the energy of the photon ) was discussed in 
detail already by Gottfried and Yennie fl4|| p|, 

2 When only one mass is excited, as it is in the present case, the phase factor depends 
only on the first and last points of the interaction. 
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The starting point in the derivation of Eq. (Q) is the generalization of 
the Mueller-Patel formula [g] for triple-pomeron coupling. In the case of N 
nucleons it reads 

F^(r,n,...r N ;Y,y) 

= e i( Zl - ZN )mx P f (n^^k r ( Xk , x ^) ni (r k ,x k ,y*)) 
J \k=l Xk X k / 

xn N (r,x' 1 ,...,x' N ;Y-y*,y-y*). (6) 

Here 

c+ioo 
c—ioo 

are the densities of the dipoles of the transverse size integrated over the 
trasverse positions of the dipoles within an onium of size r& moving with 
rapidity y*. n^ir, x\, . . . , x n ; Y — y*, y — y*) is the density of N dipoles of 
the transverse sizes x±,... xn again integrated over the transverse positions 
of the dipoles within the cascade of the transverse size r evolved from the 
incident onium. t(x, x') is the forward scattering amplitude of two dipoles 
of transverse sizes x and x' (in the two-gluon exchange approximation it is 
energy independent Q). The arguments n, . . . , r/y are the transverse sizes of 
the onia representing the target nucleons. Y, y are defined above in Eq. (||). 
In the terminology employing the pomeron concept, y is the rapidity which 
divides the process into the region of one pomeron (rapidities between Y 
and y) and N pomerons (rapidities smaller than y) in the N + 1 pomeron 
coupling. 

To obtain a formula for n^{r, x\ , . . . , x n ; Y — y*, y — y*) one has to solve 
an integro-diferential equation which one obtains by iV-fold differentiation of 
the generating functional given in ||. We do this by the method described in 
J10| ] for double-dipole scattering. An extension of this procedure to arbitrary 
N is described in Appendix A. At this point it should be emphasized that 
a solution in compact form, necessary to derive the result of Eq. (|3|) , can be 
obtained only for distributions integrated over the transverse position of the 
dipoles with respect to the original onium. The solution for arbitrary dipole 
position is much more involved. Fortunately, as already mentioned in the In- 
troduction, for large nuclei the transverse size of the incident onium is much 
smaller than the nuclear diameter and thus we only need the distributions 
integrated over the transverse positions inside the onium. 

In order to obtain the contribution to the shadow (Eq. (Q)) one has 
to add the contributions from = 2, 3, ... A nucleons and average them 
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over the positions of the nucleon inside the nucleus and over the virtual 
photon wave functions. Finally we have to integrate over the excited mass M 
(i.e. the rapidity y). Using (|2|) and 

a tot = 2ReF (8) 

one obtains 

C dh (A,x Bi ,Q 2 )= ~ 2 

Aa tot (n) 



( Y A 
cRe fdy< 

Vo N = 2 



The average (ipQ \ Fffi \ ipg) is the integral of multiplied by light-cone 
photon densities |l3| given by 

\9 Q (r,„)\* = ^d W (r, V ,Q), (10) 

W T (r, V ,Q) = ^ + (l- v ) 2 }Q 2 K 2 (Qr), (11) 
W h (r,ri,Q) = 2 V (1 - V )Q 2 K 2 (Qr), (12) 

where Q = [7/(1 — ry)] 2 Q, a em = , e 2 is the sum of the squares of the quark 
charges. 77 is the light-cone momentum fraction of one of the quarks in the 
photon. 

The relevant integrals J d 2 rdrj . . . can be found in [16| and they result in 
the following prescription: r 7 in (||) should be replaced by the expressions 

where 



TTr ^_(2-J)(l + ^ 



^(7) = K ^ 2 > ( i4 ) 



I L ( 7 ) = 2 . (15) 

Averaging over the nucleon positions reduces to multiplying (ipg \ F^ T \ tpq) 
by N single particle nucleon densities p(b, z^) and integrating over z^s kept 
ordered: 

m 1 ft 1 i/}q}) = ^ Q 1 Ff 1 i> Q ) 
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-oo +00 

x / d 2 b J dzi J d^e i ( zl - ZN ) mxp p(b,z 1 )p(b,z N ) 

—00 z\ 

x J dz 2 p(b, z 2 ) ■ ■ ■ J dz N ^ip{b,z N ^i) 

Z\ ZjV-2 

+00 +00 

= (if)Q I I Vq) / d 2 b f dz x f dz N e i ^^ m ^p(b,z 1 )p(b,z N ) 



-OO Z\ 

N-2 

X jt^ om I / dzp(b,z) I . (16) 



Z N 



(N-2)\ 



Inserting ( |l6[ ) into (^) we obtain the required formula for the shadow. 

For the sake of completeness let us also give the explicit formula for 
o"tot(n): 

o-tat(n) = a(r , Y) = 2vra 2 r 2 n cff j ^ V(7H( 7 ) , (17) 
where r 7 is given by Eq. (|l3|). 

3 Forward 7*-nucleus amplitude from the quasi- 
elastic scattering of a QCD dipole in the nucleus 

As is well-known since the seminal paper of Stodolsky [|l5[] and of Gottfried 
and Yennie fl4|| , the multiple elastic scattering of a photon fluctuation inside 
a nucleus contributes in an essential way to the nuclear shadowing. In 
the QCD dipole picture this corresponds to quasi-elastic scattering of an 
onium whose transverse size is distributed according to the light-cone wave 



functions of the photon given in [13]. In the present section we derive the 
formula for this contribution to the shadow. 

The iV-fold scattering contribution to the forward 7*-nucleus amplitude 
is built as a product of N onium-nucleon amplitudes (including the eikonal 
phases) sandwiched between the initial and final state of the virtual photon 
and summed over all intermediate states of the onium as follows 

(ipQ I F N l I ^0) = n efi(Q I T(z N ) I /cat_i) 
d 2 k]y-l d 2 kN~2 



N-l I T(z N -i) I k N - 2 > 



4n 2 47T 
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x{k N - 2 I T(z N ^2 I k N --i) . . . -j^rih I T ( z 3) I k ^~^~T 

x(k 2 \T(z 2 )\k 1 )^(k 1 \T(z 1 )\Q). (18) 

Here k\, . . . /ctv_i are the relative transverse momenta of the quark and an- 
tiquark forming the onium. To shorten the expression, we skip the integral 
symbols (and we shall continue to so doing). 

At this point it is important to realize that the amplitudes entering ( |l8| ) 
depend on the longitudinal position, acquiring the eikonal phase. So we can 
write 

(k' | T(z) | k) = e izk L(k' | T | k)e~ izkL , (19) 

where /cl and k' L are the longitudinal momenta of the onium before and after 
scattering. 
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The first step is to transform the transverse momentum amplitudes 
(k' | T | k) into transverse position amplitudes, so that we can use the 
explicit expression for the onium-onium forward elastic amplitude, T(r), 
derived in [§, |§. We shall use here the path integral representation, given 
by: 

c+ioo . . 

T(r,r ,Y)=na 2 rr J ^ {£) " e A ^ Y h(X), (20) 



where r and tq are the transverse sizes of the colliding onia. In terms of 
T(r,ro,Y) = T(r) we thus obtain 

(k' | T(z) | k) = e iz ( k2 - k ' 2 ^d 2 pe i(k '" k)p T(p) (21) 

with ^ 

* = 2^(1 - v y (22) 

In arriving at ( f22|) we have used the high energy approximation for viz. 

M 2 k 2 

k L = v = v ; -. (23) 

2u 2urj(l-rj) ( ' 

Similarly we obtain 

(k | T{z) | Q) = e iz ^ m -^d 2 p^ kf, {p [ T | Q), (24) 
where (p \ T \ Q) includes the photon wave function: 

(p\T\Q)=T(p)V Q (p). (25) 
Putting this in flT8| ) and rescaling the variables 

T TYIZ 

k = Qk, p=-, Q = — (26) 

Q vi 1 ~ V) 

we obtain 

(lP Q I Ff I $ Q ) = ^e i ( z i- z N)mx Bj ) n N ( Q | T | rN)d 2 rN 
x <P[r N ^i - r N ,x B j((N ~ (N-l)]T(r N -i)d 2 r N -i 

x <P[r N ^ 2 ~ rjv-i,»Bj(CiV-i - (N-2)]T(r N ^ 2 ) ■ ■ ■ &[r2 - ?"3,£Bj(C3 - C2)] 
x T^dV^n -r 2 ,x Bj (( 2 -Ci)]T(n)d 2 ri(ri \ T\Q), (27) 

with 



s 



where we have assumed that a = Re (a) + is to give a definite meaning to 
the integral. 

Further calculations are continued in the limit xbj <C 1. In this limit the 
following formula can be derived by the saddle point method for a,a'cl 



J d 2 r&(s - r, a)F(r)@(r - s' , a) 
= <P(s - s',a + a')F(s) exp 



where 



and 



a's + as' 
a + a 



aa 



(29) 



(30) 



(31) 



a + a' 

By repeated application of this formula we finally arrive at the following 
result 

l 

(■0Q I tff \^Q) = drj d 2 p | &Q(p,rj) \ 2 [n cS T(p, r , x B j)] N 



o 



-fV N , Z j(M?-M?, .) 

where T is given by (^) and 

M 2 = M 2 N+1 = -Q 2 , M? = (, N ~ 3 tu? 4 + V log[^(p)] ) ,j = 2, 



J \{r](l - rf)) 1 / 2 p 2 



Some details of these calculations are given in the Appendix B. 
The shadow for the quasi-elastic process is now 



(32) 

.,N. 

(33) 



C^(A,x B] ,Q 2 



Aa t (n) 



Re £(-1 



\N=2 



(A-N)l 



(34) 



with the following averaging over the nuclear densities 

l 

x / d 2 b 







Z N >Z N -l>...>Zl 



drjd 2 p | <pQ(p,r]) \ 2 [n cff T(p,r ,a;Bj)] J 
dz N dz N ^ 1 ....dz l p(b,z N )e~^ ZN{M2 v- M "+i ) 
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x P (b,z N ^)e-^-^ u ^- u ^... 

x . . . p(b, z 2 )e-^ M *- M ^p{h, Zl )e-^ Z1 ( M ?- M ^. (35) 

4 Conclusions and outlook 

We have shown that the QCD dipole picture can provide a well defined 
and compact formulae for nuclear shadowing of the virtual photons. The 
formulae are well suited for numerical evaluation. They contain the following 
parameters: 

(i) The pomeron intercept Ap\ 

(ii) The nucleon size parameter ro; 

(Hi) The effective number of dipoles in the proton n e g; 

(iv) The effective number of flavours (reflecting on e|); 

(v) The scales in the elastic and in the diffractive 7*-proton scattering. 

All these parameters can, in principle, be determined from the fit of the 
dipole picture to the proton data. Such a fit has been completed for the total 
7*-proton cross-section [§. That gave A P = .285, Q = ^ = .622 GeV; 
n e fje£ = 3.8; Y = log So, as long as the fits to the diffractive production 
are not available, we are left with two not fully determined parameters: the 
effective number of flavours and the scale in the diffractive dissociation. 

It is important to remember that the nuclear shadowing effects are de- 
termined by the forward diffractive amplitudes. Our formulae can thus be 
employed to obtain the cross section for the forward diffractive dissociation 
of the virtual photon on one nucleon. Indeed, such cross sections on one 
nucleon can be obtained from the formulae of Sections 2 and 3 for F^L 2 an d 

r N=2' 

p Q \*f£*(*i = *»)\*Q) > ( 36 ) 



and 



dycPpt lpt=0 (2tt) 

■(VQ\*%Ute = *2)\*Q), (37) 



d^ e \ 1 ,„ qe i 



d 2 Pt lpi=u (2tt) 2 

where pt is the transverse momentum of the final proton, and y defined 
in Eq. (I|) contains M, the diffractively produced mass. Therefore (^Fj) is, 
in fact, a differential cross section for production from the incident virtual 
photon of an object of mass M. On the other hand ( |37|) is the total diffractive 
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dissociation cross section integrated over diffractively excited masses. One 
can resolve (|37|) into the differential cross sections ^Tji \p t =o applying a 



straightforward procedure. This will be discussed elsewhere [19]. 

Our formulae, given in Sections 2 and 3 refer to two mechanisms de- 
scribed in the Introduction: (a) the quasi-elastic interaction of the qq dipole 
and (b) the multiple-pomeron coupling. One must remember, however, that 
these two mechanisms are not mutually exclusive: they mix in the multiple 
scattering terms, i.e. there are collisions in which they both take place. Such 
mixed amplitudes are also readily calculable by the methods developped in 
the present paper. 

Any further discussion strongly depends on the outcome of the numerical 
estimates. Before this is available let us simply list the problems which, in 
our opinion, seem interesting. 

(a) The relative importance of the multi-pomeron and the quasi-elastic 
interactions. Apart from its primary interest, it is essential for an estimation 
of the importance of the mixing terms. 

(b) An extrapolation for very small XBj where the unitarity is expected 
to break down. It is likely that for nuclear targets this effect will occur much 



earlier (i.e. for larger XBj) than estimated in 17]. 



(c) The dependence of the kind a; B j p which appear in BFKL forward 
amplitudes imply the existence of a large ratio of the real to imaginary part 



of the forward amplitudes. We have shown already |7], 18] that this influences 
significantly the predicted amount of shadowing. It will thus be interesting 
to discuss this problem again with the present, more precise, formulation. 



Appendix A 

In this Appendix we indicate some details of the derivation of the formula 
for nTv, Eq. (|). 

We start with n<i- 

Using the generating function given in || and assuming that evolution of 
the two exchanged pomerons goes until the rapidity y is reached, we obtain 
the following equation for ri2 (c.f. also Eq. (52) of |2|): 
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dn 2 (x m ,Y,y,x,x') 



dY 

bi ' 



aNc_ I x ±d ^ ni ( Xl2j y jX } ni ( X02j y jX >}s(Y - y) 



7T J Xf\()X 
R 



02^12 



+ 2aNc I dxi 2 K(x i,x 12 )n 2 (x 12 ,Y,y,x,x') , (38) 



7T 

where n\ is the single dipole density and 

J_ r x% 1 d 2 x 2 
2vr J 



K(x i,xi 2 ) = — I — — 2 ^( Xo1 ~ Xl2 ) io s(xoi/p) (39) 

T 02 X 12 



is the Lipatov kernel (c.f. |j|]) which satisfies the eigenvalue equation: 

K(x 01 ,x 12 )xi 2 dx 12 = SqixM ( 40 ) 



with x(A) given by (H). We now introduce Mellin transforms to take advan- 
tage of ©: 

c+ioo 

n 2 (^oi, Y,y,x,x') = J -^-:h 2 ('y,Y,y,x,x')xl 1 (41) 

c—ioo 
oo 

n 2 (j,Y,y,x,x') = / dx^x^" 1 n 2 (x i, Y, y, a;, a;') . (42) 



Substituting (|4lj) into (p8|) and using ([40|) we obtain a differential equation 
for 77,2 which can be easily solved giving 

h 2 (~f,Y,y,x,x') =e^( Y -^g 2 ( 7 ,y,x,x'), (43) 

where ^(7) is given by (0) and 



2 ,2 



, aiV c /• _i_ 7 /■ sc^tTa^ . « , 

92{l,y,x,x) = —2- / dx ix 01 / — 5 — — 77i(xi 2 ,y,x)ni(x 2,y,x )• 

/? uz ^ 

(44) 

Substituting (g|) and (||) into (§]]) we finally have 

c+ioo 

n 2 (r,Y,y,x,x') = ^ f ^l r 7 e ^)(Y-y) 

7T Z 7 27T7 



C— lOO 



x / dxoix^ 1 7 / X ° 2 ld 2 ^- ni(x 12 ,y,x)n 1 (xo2,y,x'). (45) 

x 02 x 12 
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We now turn to TI3. 

Using the generating function given in || and assuming that evolution 
of the three exchanged pomerons goes until the rapidity y is reached, we 
obtain the following equation for n^: 

dn 3 (x 01 ,Y,y,x,x',x") 
dY 

= — T I X °2 T \ n i( x i2,y, x)n 2 {xQ 2 , Y, y, x' ,x")] sym 5(Y - y) 

+ 2aNc [ dx 12 K(xo 1 ,x 12 )n 3 (x 1 2,Y,y,x,x',x"), (46) 
vr J 

where n2{xQ2,Y,y, x',x") is given by (|45|) . The symbol sym denotes a sum 
of three terms, symmetrized with respect to X . X • X • 

Introducing Mellin transforms, as in (f|l|), (|42|), substituting them into 
(E§) and using ( |40D we obtain a differential equation for h 3 . The solution is 
found again in the form 

7i 3 (7, Y, V, x, x', x") = e^W( Y -y)g 3 ( 7 , y, x, x', x"), (47) 

where 



aN r 



93{l, y,x,x',x") = J dx ix Q l 7 



f x^ d^x 2 

x / 2 x 2 [n 1 (xi2,y,x)n2(x 2,y,y,x',x")] sjm , (48) 



Xnn X 



02-^12 



and n 2 (x 2,y,y,x',x") is given by (||). 

The formulae for n.4, 77,5, . . . are derived in a similar way. 

The integrals over x\, . . . xjy; x[, . . . , x' N in (g) are performed using the 
formula for the forward onium-onium amplitude T(r,r',y) 

f dx dx' 

T(r,r',y) = / -r(x, x )ni(r, y*,x)ni(r', y - y*, x') 

J x x 

= vraV J ^ A[l)Y Qj) ^Mt) (49) 

c— ioo 

which was first derived in |§ and g, |] in a sligthly more simplified form. 
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Using (E9) and (H5J), (H) gives 



c+ioo 



where 



R^n^^y) = e^i-^W^ ri r 2 f ^ e A M (Y-y) 

TT Z J 2tT\ 

c— ioo 

c+ioo 

c— ioo 
c+ioo 

/ f^-V^H^^A^), (50) 



^(7,Ai,A 2 )= / dx ix 01 1-7 / °2 2 - aQ2~ Al a;i2~ A2 , (51) 

x 02 x 12 



J2(7, Ai, A2) can be calculated using 



z 01 d x 2 = dxwdxi2 f kdkj Q (k XQ1 }j Q (kx 02 )Jo(kx 12 ), (52) 

xyx\ 2 x 02 X12 J 



from which the following useful identity can be obtained: 

00 00 

f r -i- 7 f x 2 01 d 2 x 2 \ v G(A)G(A') o 

/ dxoix 01 / 2 2 ^02^12 = ^ — TTT^ — 2vrz5(7 - A - A ). (53) 

J J X Q2 X 12 ^wJ 



R 



Using (|53|) in (pQ) and setting ri = r2 = ro we obtain the formula (||) 
for -F^ 11 - With the same technique one can obtain (Q) for more than two 
collisions. We indicate below the derivation for three collisions. 
Substituting (HEf) into (S3) and using fl47|) we obtain 



2 c+ioo 



n 3 (x 01 ,Y,y,x,x',x") = J 7^.xli<i A ^ ){Y y) J dx ix 01 

c— ioo 
2 £, c+ioo ^ t 



r> t 02 J '12 
it c— ioo 

2 2 

„' (J \-l-7' A ( X 0l) ^ X 2 



x / ^01(^01) 1 7 / / / 01 x 2 /^ l 2 M x 'i2, V, x')ni(x' 02 , y, x")] sym . (54) 

J J \ X 02J \ X 12) 
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When (^4|) is substituted into (y) with n = r 2 = r% = r$ we obtain 

Ft(r,r ;Y,y)=e i ^-^ m ^ 

c+ioo oo 



V 7T 



x3' * ' 



9 ,9 C+i °° t / 



X 02 X 12 J Z7TZ 

K c— loo 

/ A 2. 2 / 



x f dx' m (x' m )- 1 ^' f ^^ T(x' 12 ,r ,y)T(x' 02 ,r ,y), (55) 

{ J \ X Q2) \ X 12) 

U H, 

where T(r, r',y) is the dipole-dipole forward scattering amplitude given 
by @. 

The integrals over dxo±,d 2 X2 and dx' 01 d 2 x 2 can be performed using again 
the identity (|5^) and taking into account the fact that the dependence of 
T(r, r', y) on r and r' is in the form of a power law. The result of these 
operations is the Eq. (||) for N = 3. 



Appendix B 

In this Appendix we give a few intermediate steps of the calculations of 
Section 3 which lead to Eq. (|32|). 



First, by repeated application of (29) we obtain 



<V>Q I Ff I i, Q ) = n&Q-Vfr-^^Q I T I r N )d 2 r N 

x ^(n - rAr,rr B j(CAr - CAf-i))5 , (r A r,ri;CAr, • • . ,(i)d 2 ri{n \ T | Q),(56) 



where 



5(r JV ,r 1 ;Civ,...,Ci)= II 2^( A i' r °' F ) r i' 

i=2 

. jv-l 

x ex P I + i2 2^ X 3 X k±3^jk ) (57) 

3 j,k=2 3 k 



with 



(z N - Zj)fi + (zj - zi)r N 
z n - z\ 



(58) 
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(z N - Zj)(Zj - Zj) 

Zn-Zl 
(z N - Zj)(z k - Zl) 



Zn ~ Zl 

A(X,r ,Y) = n cS 7ra 2 r 2 - x e A ^ Y h(X) 



(59) 

(60) 
(61) 



is taken from the integrand of Eq. (|20| ) and n c g added for the reasons ex- 
plained in the main text. 

We finally perform the integral over d 2 r^, again in the saddle point 
approximation which is valid for small xbj- Then <&{r\ — rjy, xbj(Ctv — Ci)) is 
very strongly peaked around r n ~ r±, so that in leading order all r^-'s become 
r\. This approximation puts a new kind of factor, Vipq/ipQ, into the over 
all phase of the expression and we obtain (returning to Zj's and pj's) 



I Ff | Vq) = nf ff e ix Bi(^-N) / drj j d 2 p j ,2 
1 



x exp 



N 



2^ 



N 



rj(l - rj)p 



3=2 



V 



+ (z n - Zl, 



^q(p) 2 



3=1 



2tti 



A{\ v r ,Y)p x i . 



(62) 



One can check by a direct substitution of M 2 given by (|33|) into (32) that it 

reduces to (|62]) . The version of (tpQ | F^ cl | iJjq) given in (^3|) has, however, 
a clearer physical interpretation. 
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